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Abstract. New observations of eleven compact galaxies in the 12 CO J = 2 — 1 and J = 3 — 2 transitions 
are presented. From these observations and literature data accurate line ratios in matched beams have been 
constructed, allowing the modelling of physical parameters. Matching a single gas component to observed line 
ratios tends to produce physically unrealistic results, and is often not possible at all. Much better results are 
obtained by modelling two distinct gas components. In most observed galaxies, the molecular gas is warm (Tk = 
50 — 150 K) and at least partially dense (n(H2) > 3000 cm -3 ). Most of the gas-phase carbon in these galaxies 
is in atomic form; only a small fraction (~ 5%) is in carbon monoxide. Beam-averaged CO column densities 
are low (of the order of 10 16 cm -2 ). However, molecular hydrogen column densities are high (of the order of 
10 cm - ) confirming large CO-to- H2 conversion factors (typically X = 10 21 - 10 22 cm -2 / K kms -1 ) found 
for low-metallicity environments by other methods. From CO spectroscopy, three different types of molecular 
environment may be distinguished in compact galaxies. Type I (high rotational and isotopic ratios) corresponds 
to hot and dense molecular clouds dominated by star-forming regions. Type II has lower ratios, similar to the mean 
found for infrared-luminous galaxies in general, and corresponds to environments engaged in, but not dominated 
by, star-forming activity. Type III, characterized by low 12 CO (2-l)/(l-0) ratios, corresponds to mostly inactive 
environments of relatively low density. 
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1. Introduction 

Compared to large spiral galaxies, dwarf and compact 
galaxies are difficult to detect in molecular lines, even in 
those of the relatively abundant CO molecule. However, 
molecular clouds constitute the unique environment out of 
which the stars are formed and knowledge of their occur- 
rence, spatial distribution, mass and physical condition is 
essential to gain insight in the process of star formation in 
galaxies. For this reason, many surveys of dwarf and com- 
pact galaxies have been conducted over the last decades, 
although generally with low detection rates (Tabled . 

Detection rates are highest for 'tidal' dwarfs (cf. Braine 
et al. 2001), 'big' dwarfs (Sm) and galaxies that are only 
modestly metal-poor. They are lowest for 'small' dwarfs 
(Im), blue compact dwarf galaxies (BCDGs) and practi- 
cally non-existent for very metal-poor dwarfs (Taylor et 
al. 1998; Barone et al. 2000). In the extensive survey by 
Albrecht et al. (2004), detection rates are 72% and 39% 
for Sm and Im galaxies respectively (see also Leroy et al., 
2005). 

Velocity-integrated 12 CO(l — 0) fluxes are often com- 
bined with the so-called CO-to-H2 conversion factor X to 
derive a beam-averaged H2 column-density, hence an es- 



Table 1. CO surveys of dwarf and compact galaxies 



Reference 



Number 
Obs. Det. 



Line 



Elmegreen et al. (1980) 


7 


1 


1- 





Gordon et al. (1982) 


7 





1- 





Israel & Burton (1986) 


12 





1- 





Thronson & Bally (1987) 


22 


10 


1- 





Tacconi & Young (1987) 


15 


6 


1- 





Arnault et al. (1988) 


12 





1- 





Sage et al. (1992) 


15 


8 


1-0/2- 


1 


Hunter & Sage (1993) 


5 





1-0/2- 


1 


Israel et al. (1995) 


25 


6 


1- 





Gondhalekar et al. (1998) 


29 


5 


1- 





Taylor et al. (1998) 


11 


3 


1- 





Barone et al. (2000) 


10 


2 


1-0/2- 


1 


Meier et al. (2001) 


8 


8 


3- 


2 


Braine et al. (2001) 


10 


8 


1-0/2- 


1 


Boselli et al. (2002) 


6 


1 


1- 





Albrecht et al. (2004) 


64 


41 


1-0/2- 


1 


Leroy et al. (2005) 


121 


47 


1- 
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timate for the H2 mass in the beam. It is widely (but not 
yet universally) accepted that this factor AT is a function 



2 



F.P. Israel: Molecular gas in compact galaxies 



of environmental conditions (such as metallicity, radiation 
density etc.). However, the extent to which X is sensitive 
to e.g. changes in metallicity is still a matter of debate. 
Arguments for a strong dependence of X on metallicity as 
traced by the oxygen abundance [0]/[H] have been pre- 
sented by Israel (1997; 2000), Barone et al. (2000) and 
Boselli et al. (2002). 

Obviously, it is much preferable to determine actual 
molecular column densities and masses from an analysis 
of the physical condition of the gas than from an empiri- 
cally determined factor that relates column density to the 
flux of an optically thick line. In principle, multi-line ob- 
servations can be used as a diagnostic tool to estimate the 
physical condition of the gas, but this requires that the 
various transitions have been observed at the same celes- 
tial position with matched beamsizes. A glance at Tabled 
shows that only a few surveys sample CO line transitions 
other than J=l-0. In a number of surveys, J=2-l and 
,7=1-0 measurements were made simultaneously with the 
same telescope, therefore with beam areas differing by a 
factor of four. Mapping the J=2-l line over the extent of 
the J=l-0 beam would make it possible to compare the 
two transitions, but this is rarely done. Only the very re- 
cent survey by Albrecht et al. combines measurements in 
the two transitions obtained with the IRAM 30m and the 
SEST 15m apertures, providing matched (SEST) ,7=2 1 
and (IRAM) J=l-0 measurements. In all preceding sur- 
veys, including the CSO 12 CO (3 — 2) survey by Meier et 
al. 2001, line ratios could only be estimated as a function 
of (assumed) source extent (point source, extended source, 
or inbetween). As can easily be seen, this introduces very 
large uncertainties, effectively ruling out meaningful mod- 
elling of physical paramters. 

In this paper we present new observations, which to- 
gether with measurements published in the literature, al- 
low us for the first time to construct reliable line ratios in 
matched beams for a sample of about a dozen dwarf or 
compact galaxies. 

2. Observations 

All observations described in this paper were made with 
the 15m James Clerk Maxwell Telescope (JCMT) on 
Mauna Kea (Hawaii) 1 , mostly between July 1995 and 
December 1996. They were made in beamswitching mode 
with a throw of 3' in azimuth using the DAS digital auto- 
correlator system. When sufficient free baseline was avail- 
able (i.e when the detected line was sufficiently narrow), 
we subtracted second or third order baselines from the 
profiles. In all other cases, linear baseline corrections were 
applied. All spectra were scaled to a main-beam brightness 
temperature, T m b = T^/rj^. We used values of r) m b appro- 

1 The James Clerk Maxwell Telescope is operated on a 
joint basis between the United Kingdom Particle Physics and 
Astrophysics Council (PPARC), the Netherlands Organisation 
for Scientific Research (NWO) and the National Research 
Council of Canada (NRC). 



priate to the epoch of observation. These values were be- 
tween 0.64 and 0.72 at 230 GHz, between 0.53 and 0.60 at 
345 GHz, and between 0.50 and 0.53 at 461 GHz. Spectra 
of observed positions are shown in Fig.^ We also made 
small maps of the four galaxies (IC 10, Haro 2, NGC 4194 
and NGC 6052) distinguished by relatively strong J=3-2 
CO line emission. These maps are shown in Figs.El&El 

3. Results 

3.1. Detections 

Details of the observed galaxies and the observational re- 
sults obtained are summarized in Tabled In this table, 
Column 1 gives the names, with alternative names in 
the notes at bottom. Right ascensions and declinations in 
columns 2 and 3 are those of the observed (0,0) position. 
The receivers were tuned to the radial velocity given in 
column 4; as is obvious from Fig. ^ this was usually very 
close to the CO velocity. Column 6 gives the integrated 
absolute blue magnitude as retrieved from photometric 
data in the NASA-IPAC NED and calculated for the as- 
sumed galaxy distance listed in Column 5. there. Columns 
7 and 8 identify the observed 12 CO transition and the 
concomitant FWHM beamsize. Finally, columns 9 and 10 
present the observed peak main-beam brightness temper- 
ature, and the CO intensity integrated over the observed 
profile. Quoted errrors are r.m.s. values; upper limits are 
three times the r.m.s. value. 

As Fig.n and Tabled show, emission in the J=2-l 
12 CO transition was detected from all observed galaxies, 
although the emission from 2 Zw 40 was close to the de- 
tection threshold. 2 Zw 40 was not detected in the J=3-2 
transition nor was NGC 2537 in spite of its clear detec- 
tion in the J=2-l transition. By convolving the maps of 
IC 10, Haro 2, NGC 4194 and NGC 6052 to a resolution of 
21", we also obtained for these galaxies the J=3-2 12 CO 
intensity matched to the J=2-l beamsize (Tabled- 

3.2. Line ratios 

We have determined rotational line ratios by directly com- 
paring transitions observed with beams of identical or very 
similar size, using our own measurements and results from 
the literature, summarized in Table|3| Where observations 
in identical beams were lacking, we extrapolated CO fluxes 
to matching beamsizes. This was done either directly by 
plotting observed CO flux (TablesEl & EJ) as a function 
of beamsize and reading off the desired flux value at the 
relevant beamsize, or indirectly by using the equivalent 
CO source size. The equivalent size can be determined 
from observations with different beams (and not neces- 
sarily in the same transition) under the assumption that 
the CO emission has a radial gaussian brightness distri- 
bution with circular symmetry. It can also be estimated 
from high-resolution (array) maps of the CO distribu- 
tion. Actual methods used are found in Sect. 3.3 dealing 
with the individual galaxies and the resulting ratios are 
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Fig. 1. J=2-l and J=3-2 12 CO emission-line spectra observed towards the sample dwarf and peculiar galaxies. 
Horizontal scale is Vlsr in kms -1 , vertical scale is main-beam brightness temperature T m b in Kelvins. 



given in Table^J together with metallicities [0]/[H] taken ies observed by us, whereas the remainder consists of data 

from data compiled by Taylor et al. (1988) and Meier et culled from the literature. The bottom section contains 

al. (2001) and the IRAS far-infrared 60/xm/lOO/im flux- galaxies experiencing starburst activity, 
density ratios taken from Lonsdale et al (1989). Tabled Not included in Tabled is the dwarf galaxy UM 465 

consists of three parts: the top section contains the galax- for which CO measurements imply a (2-l)/(l-0) ratio of 
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Table 2. JCMT 12 CO observations of compact galaxies 
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Note: (a) Mkn 86 = Arp 6; (b) Mkn 33 = Arp 233; (c) Mkn 201 = 1 Zw 33 = Arp 160; (d) 1 Zw 89; (e) Mkn 297 = Arp 209. 




Fig. 2. Contour maps of integrated CO emission. Lowest contour is always at zero. Left: IC 10 ( J=3-2) with main- 
beam brightness temperature contours in multiples of 1 K. Center: NGC 4194 (J=2-l) with contours in multiples of 
4 K. Right: Haro 2 (J=3-2) with contours in multiples of 0.75 K. 



1.16 ± 0.29 (Sage et al. 1992; Barone et al. 2000) and re- 
cently published results by Albrecht et al. (2004) which 
allow direct determination of the 2 — 1/1 — ratio for an- 
other six dwarf galaxies: NGC 145 (0.74+0.13); NGC 2730 
(0.62+0.17); NGC 4532 (0.59 ± 0.12); NGC 6570 (0.41 ± 
0.12); NGC 7732 (0.49 + 0.24); IC 3521 (0.49 + 0.14) as 
well as upper limits for four more: NGC 178 (< 0.77); 



NGC 1140 (< 1.36); NGC 3659 (< 0.17); NGC 4234 
(< 0.80). 

On average, CO emission from compact galaxies is 
weak. As emission from the 13 CO isotope is usually at 
least an order of magnitude less than that from 12 CO, very 
few galaxies have been measured in this line. In Table 0we 
list the available information on the galaxies from Table0J 
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Table 3. Summary of literature 12 C0 data on observed galaxies 
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References: 1. Israel, Tacconi & Baas (1995); 2. Tacconi & Young (1987); 3. Thronson & Bally (1987); 4. Gondhalekar et al. 
(1998); 5. Albrecht et al. (2004); 6. Sage et al. (1992); 7. Meier et al. (2001); 8. Barone et al. (2000); 9. Sofue et al. (1990); 10. 
Yao et al. (2003); 11. Aalto & Hiittemeister (2000); 12. Casoli, Dupraz & Combes (1992); 13. Devereux et al. (1994); 14. Taylor, 
Kobulnicky & Skilman (1998); 15. Thronson et al. (1988). 

Table 4. Integrated rotational line ratios in compact galaxies 
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In effect, this sample is limited to actively star-forming 3.3. Individual galaxies 
galaxies. 

NGC 55 is a small SB(s)m galaxy, roughly as we would 
expect to see the LMC edge-on, member of the nearby 
Sculptor Group and rich in interstellar gas. Dettmar & 
Heithausen (1989) used the SEST to make a small J=l-0 
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Fig. 3. Grid map of the J— 3-2 CO emission from 
NGC 6052. Horizontal scale is radial velocity Vlsr in 
kms , vertical scale is main-bean brightness tempera- 
ture T m b in Kelvins. 



12 CO map of the brightest part. At the position of the 
CO peak in this map, Becker & Freundling (1991) also 
measured the J=l-0 13 CO intensity. No further CO data 
on this interesting galaxy have appeared in the literature. 
Our J— 2—1 spectrum (not shown) was taken very close to 
the CO peak, at the same position as our earlier J=l— 
spectrum (Israel, Tacconi & Baas, 1995). 

IC 10 (dlrr IV) is an unusual Local Group galaxy; 
it is one of the most luminous dwarf galaxies in the far- 
infrared (Melisse & Israel 1994), and appears to be the 
nearest representative of the class commonly referred to 
as Blue Compact Galaxies (BCG). It contains two major 
groupings of HII regions (see Hodge & Lee 1990; Yang 
& Skillman 1993; Chyzy et al. 2003), the brightest of 
which is alternatingly known in the literature as IC 10A 
or IC 10-SE. The map in FigJ2 shows the CO distribu- 
tion associated with IC 10-SE, which resembles that of 
far-infrared and submillimeter emission from heated dust 
(Thronson et al. 1990) and that of ionized carbon ([CII]) 
emission (Madden et al. 1997). Studies of the dense molec- 
ular medium in IC 10-SE were published by Petitpas & 
Wilson (1998) and Bolatto et al. (2000). The (2-l)/(l-0) 
ratio in Tabled] was taken from the former; the (3-2)/(2- 
1) ratio is the mean of their determination and our higher 
value. 

NGC 1569 (IBm, Sbrst) is an extreme starburst 
dwarf galaxy (Israel 1988). Weak CO detected by vari- 
ous authors (Greve et al. 1996; Taylor et al. 1998, 1999; 
Meier et al. 2001; Albrecht et al. 2004) has been mapped 
and discussed in Muhle's 2003 Ph.D. thesis from which 



we took the ratios in Tabled] (average of values from her 
Table 3.9). The extent of CO emission is slightly less than 
30". 

LGS 3 (dlrr/dSph), a Local Group galaxy in Pisces, 
is one of the faintest dwarf galaxies known. It was sur- 
prisingly detected in CO by Tacconi & Young (1987) with 
a 45" beam. Our measurements were made at the same 
position, but with smaller beams. The range of ratios in 
Tabled was obtained by assuming the CO source to be 
pointlike or extended respectively. 

2 Zw 40 is a gas-rich galaxy alternately classified as a 
BCG or an HII merger galaxy (Beck et al. 2002, and ref- 
erences therein). CO was detected by Thronson & Bally 
(1987), Tacconi & Young (1987), and Sage et al. (1992). 
Our spectra were taken at the position used by the lat- 
ter. The integrated intensities as a function of beamsize 
suggest an equivalent CO size of 23" which was used to 
obtain the (2-l)/(l-0) ratio in TableU 

He 2-10 (IOpec, Sbrst) is a dwarf galaxy experi- 
encing a strong starburst caused by an ongoing merger. 
The equivalent CO size is about 12", consistent with the 
OVRO J=l-0 CO map published by Kobulnicky et al. 
(1995). The CO ratios in Table^Jwere taken from data in 
that paper, from Meier et al. (2001) and from Baas, Israel 
& Koornneef (1994). 

NGC 2537 (SB(s)m pec) is also considered to be 
a BCG and often referred to as Mkn 86 or, after its 
optical appearance, as the Bear Claw. Not detected by 
Thronson & Bally (1987), 12 CO emission was succesfully 
measured in the J=l-0 and J=2-l transitions by Sage 
et al. (1992). Our measurements were made at the same 
position, and the line ratios follow directly from the data 
listed in Tablcs|21 & 01 NGC 2537 has also been mapped 
in J=l-0 and J=2-l 12 CO by Gil de Paz et al. (2002). 
Our 12 CO (2-l)/(l-0) ratio agrees with their less certain 
value 1.06±0.40 

NGC 2976 (SAc pec, HII) is part of the M 81 group. 
Optically, it is similar to NGC 2537 and NGC 1569 in 
having a complex central structure surrounded by an 
extended and relatively featureless distribution of stars. 
Bright HII regions occur at either end of the major axis. 
CO was clearly detected by Thronson & Bally (1987) in 
a large beam (100"). Our CO measurements refer to the 
center and do not include either of the major HII region 
complexes. This central part has also been mapped in 
J=l-0 12 CO at high resolution by Simon et al. (2003) 
using BIMA. Their map shows extended emission consis- 
tent with the equivalent CO size of about 1' suggested by 
the single dish measurements. The values in Table0]were 
calculated from the data in Tables|3 & [3] by assuming a 
CO source size of 60". 

NGC 3077 (10 pec, HII) is also a small member of the 
M 81 group. CO emission from this galaxy was discovered 
by Becker et al. (1989). Our measurements coincide with 
the peak in their map. High-resolution OVRO maps of 
the galaxy in the J = l-0 and J=2-l CO transitions were 
published by Meier et al. (2001) and Walter et al. (2002) 
and show significant substructure. The ratios in TablcQ] 
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were derived from data by Becker et al. (1989), Meier et 
al. (2001) and Albrecht et al. (2004) for an equivalent CO 
source size of 25". 

Haro 2 (Im pec, HII) is another BCG also known as 
Mkn 33. CO was detected by Thronson & Bally (1987), 
Sage et al. (1992), Israel, Tacconi & Baas (1995), Barone 
et al. (2000), and Meier et al. (2001) including measure- 
ment of the J=2-l and J=3-2 12 CO transitions. An 
OVRO map of the J=l-0 12 CO distribution was pre- 
sented by Bravo- Alfaro et al. (2004) and shows that most 
of the emission originates in a source of size 15"-20" (cf. 
Fig.[2). Our (2-l)/(l-0) 12 CO ratio is based on the mean 
of the J=l-0 CO fluxes in Tabled] The implied (3-2)/(2- 
1) ratio appears to be dependent on the beamsize used. 
The mean ratio for the 21/22" beam is 0.63; whereas the 
the mean for the 12/14" beam is 1.08. In Tabled we list 
the average of these two. 

NGC 3353 = Haro 3 = Mkn 35 (Irr, HII) is a 
BCG detected by Thronson & Bally (1987), Tacconi & 
Young (1987), Sage et al. (1992), and Meier et al. (2001). 
The available data only allow a reliable determination of 
the (3-2)/(l-0) ratio, as the J =2-1 measurement by Sage 
it al. (1992) appears to be too low. 

NGC 4194 (IBm pec, SB0 pec, HII) is also known 
as the BCG 1 Zw 33, Mkn 201 or the Medusa Merger. 
Measurements of the lower three 12 CO transitions have 
been published by Casoli et al. (1992), Devereux et al. 
(1994) and Aalto et al. (2001). An OVRO J=l-0 12 CO 
map was presented by Aalto & Hiittemeister (2000). The 
equivalent CO size of 8" suggested by the single-dish 
measurements is consistent with the extent of the major 
CO concentration in their map. The data summarized in 
Tableware, unfortunately, not entirely consistent. In par- 
ticular, we have decided to ignore the high J=2-l 12 CO 
value from Casoli et al. (1992), which is inconsistent with 
all other measurements. To obtain the J=2-l/J=l-0 ra- 
tio, we have divided our J=2-l value from Tabled by 
the mean of the 22/24" J=l-0 data in TableEJ For the 
J=3-2/J=2-l ratio we have taken the mean of our own 
(21") values in TableU (0.78) and the 11/14" values in 
that Table and in Tabled (0.86). In view of the uncertain- 
ties, we have not attempted to correct the latter for finite 
source and beam sizes. 

NGC 4214 (IAB(s)m, HII) was measured in 12 CO by 
Tacconi & Young (1985), Thronson et al. (1988), Ohta et 
al. (1993), Becker et al. (1995), and Taylor et al. (1998). 
An OVRO map by Walter et al. (2000) shows three CO 
concentrations along the major axis separated by 40", each 
about 10-15" in size. Our measurements (cf. FigJ^J sample 
the southern maximum. The range of ratios in Tablc^] 
refer to source sizes varying from great extent (high ratio) 
to essentially unresolved (low ratio). The (3-2)/(2-l) ratio 
cannot be determined accurately, but the upper limit is 
physically significant. 

NGC 5253 (Impec, HII, Sbrst), also known as 
Haro 10, is a small but pronounced star burst member of 
the M 83 group. The J=l-0 observations in various beams 
(Wiklind & Henkel 1989; Turner et al. 1997; Taylor et al. 



Table 5. Integrated isotopic line ratios in compact galax- 
ies 



Galaxy 


7 i n 
J = 1-1) 

(n) 


12 l 13 r^r\ 

CO/ CO 
J =2—1 


J=3-2 


Kef 


NGC 55 


15±4 






1 


IC 10-SE 


9±1 


13±3 


12±2 


2 


He 2-10 


20±4 


20±4 


16±5 


3 


NGC 1569 


23±8 


36±9 




4 


NGC 4194 


19±4 






5 


NGC 6822-HV 


23±7 




11±3 


6 



References: 1. Becker & Freundling (1991); 2. Petitpas & 
Wilson (1998); Bolatto et al. (2000); redetermined from 
archive data; 3. Baas et al. (1994); Kobulnicky et al. (1995); 
4. Miihle (2003); 5. Aalto et al. (2001); 6. Israel et al. (2003), 
redetermined. 



1998) suggest an equivalent CO size of 35". This, together 
with the J=2-l and J=3-2 obervations by Meier et al. 
(2001) implies the ratios given in TablcQ] 

NGC 5633 ((R)SA(rs)b) was detected in CO by 
Thronson & Bally (1987). No other CO measurements 
have appeared in the literature, so that we can only con- 
strain possible ratios by the assumption of either pointlike 
or extended CO emission. Again, only the (3-2)/(2-l) ra- 
tio is sufficiently constrained to have physical meaning. 

NGC 6052 is an otherwise unclassified BCG more 
commonly known as Mkn 297 with an optical size of about 
20", possibly the result of a merger. Comparison of all 
available measurements suggests that Sage et al. (1992) 
overestimated the J=l— intensity; the ratio derived from 
14" beam observations (assuming identical source struc- 
ture in all CO transitions) is more plausible and given in 
Tabled Our J=3-2 map (Fig.|3|) shows the CO source to 
be resolved whereas the J=l-0 map by Sofue et al. (1990) 
only shows marginal resolution. Taken together, these re- 
sults suggest a CO extent of about 15". 

NGC6822 is a Local Group galaxy similar to the 
LMC. The entries in Tabled were taken from Israel et 
al. (2003) and refer to individual pointings on the CO 
cloud complex associated with the star-formation region 
Hubble V. 

4. Analysis and discussion 

4.1. Physical condition of the gas 

The observed 12 CO and 13 CO transitions can be an- 
alyzed to provide constraints on the physical condition 
of the molecular gas in the compact galaxies concerned. 
To this purpose we have used the large-velocity gradi- 
ent (LVG) radiative transfer models described by Jansen 
(1995) and Jansen et al. (1994). They provide model line 
intensities as a function of three input parameters: gas ki- 
netic temperature Tk, molecular hydrogen density n(H2) 
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Table 6. Physical parameters of model clouds 





Kin. 

Temp. 
T k 
(K) 


Gas 
Density 
n(H 2 ) 
(cm" 3 ) 


Column 
Density 
N(CO)/dV 
( cm" 2 / kins -1 ) 


Contribution 
to J=2-l 
Emission R21 


Observed 
Ratios 
R32 n 


7-2 


r-3 


R21 


Model 
Ratios 
R32 ri 


r-> 


f3 








NGC 2537, NGC 2976, Haro 2, NGC 4194, NGC 6052 
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^000 
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1.00 0.91 
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IC 10-SE 


















1 

2 


30 
100 


10 4 - 10 s 
100 


1.0 x 10 17 
1.0 x 10 17 


0.15 0.64 
0.85 

NGC 1569 


0.98 


9 


13 


12 


0.76 


0.93 11 


11 


13 


1 
2 


100 
100 


10 5 
1000 


1.0 x 10 17 
0.3 x 10 17 


0.25 1.15 
0.75 


1.05 


23 


36 




1.23 


0.99 23 


20 


18 










He 2-10, NGC 3077 
















1 
2 


60 
100±50 


10 5 
500-1000 


0.6 x 10 17 
0.3 x 10 17 


0.20 0.90 
0.80 


1.14 


20 


20 


16 


0.99 


0.83 21 


18 


17 










NGC6822-HV, NGC 5253 
















1 
2 


30 
30-150 


10 5 
10 5 


0.6 x 10 17 
0.3 x 10 17 


0.75 1.35 
0.25 


1.03 


23 




11 


1.45 


0.99 22 


14 


11 



and CO column density per unit velocity (jV(CO)/dV). 
Comparison of model to observed line ratios allows iden- 
tification of the physical parameters best describing the 
actual conditions in the observed source. Beam-averaged 
properties are determined by comparing observed and 
model line intensities. Tables^] and show that in sev- 
eral galaxies the input parameters are not fully con- 
strained by the observations. At least some constraint on 
gas conditions can be derived for NGC 2976, NGC 3077, 
Haro 2 (Mkn 33), NGC 4194 (Mkn 201), NGC 5253, and 
NGC 6052 (Mkn 297). Sufficient information for a proper 
analysis is provided only by (the star-forming complexes 
in) IC 10-SE, NGC 1569, NGC 6822-HV and He 2-10. 

4.1.1. Limitations of single-component modelling 

First, we have attempted to fit the observed line ratios 
with a molecular gas at a single temperature, a single 
density and a single velocity gradient by searching a grid 
of model intensity ratios corresponding to temperatures 
T k = 10 - 250 K, densities rc(H 2 ) = 10 2 - 10 5 cm~ 3 , and 
gradients N{CO) /dV = 6 x 10 15 - 3 x 10 18 cur 2 / km s" 1 ) 
for values matching the observed intensity ratios. Among 
the galaxies listed in TableU five (NGC 2537, NGC 2976, 
Haro 2, NGC 4194 and NGC 6052) have almost identi- 
cal line ratios with mean values (2-l)/(l-0) = 0.90 and 
(3-2)/(2-l) = 0.78. These ratios are accurately repro- 
duced (TableEJ) by a warm and moderately dense molec- 
ular gas of temperature T^ in — 60(+40, —20) K, den- 
sity ra(H 2 ) = 3000 cm- 3 , and gradient N(CO)/dV = 
0.6-1.0xl0 17 cm- 2 (kms- 1 )- 1 . The J=l-0 12 CO/ 13 CO 



ratio in NGC 4194 (TableEJ agrees well with this. The 
intensities imply a filling factor 3.3 x lO -3 and beam- 
averaged CO column-densities range from « 1 x 10 16 cm" 2 
(NGC 2976, Haro 2) to 5 x 10 16 cm~ 2 (NGC 6052) to 
1 x 10 17 cm" 2 (NGC 4194). The relatively low (2-l)/(l-0) 
ratios (mean: 0.56) for the Albrecht et al. (2004) galaxies 
(Sect. 3.2), and for 2 Zw 40 (Table|U), suggest the domi- 
nating presence of fairly low-density (n(H2) < 800 cm -3 , 
iV(CO)/dV = 0.3 - 0.6 x 10 17 cm- 2 (kms- 1 )- 1 r gas at 
undetermined temperatures. 

The physical meaning of these results should be es- 
tablished by further observations, especially of the 13 CO 
isotope, but we suspect that additional information will 
only serve to rule out single-component fits after all. It 
is perhaps telling that no other galaxy in Tabled is sat- 
isfactorily fitted by a single component. Moreover, even 
in cases where a single component appears to provide a 
good (LVG) fit to the observations, caution should be ex- 
ercised in accepting the result as a physical reality, espe- 
cially when the fit is based on limited data. This may be 
illustrated by the case of IC 10-SE. The relative intensi- 
ties of the J=l-0, J=2-l and J=3-2 12 CO lines together 
with the J=l-0 and J=2-l 12 CO/ 13 CO isotopic ratios 
are perfectly fitted by a single hot (Tk; n = 100 K) and 
tenuous (n(H 2 ) = 100cm- 3 ) molecular gas component. 
Note that this result is based on a total of four line ra- 
tios. As many determinations in the literature are based 
on three or only two ratios, this is more than the num- 
ber commonly used in such analyses. Yet, this apparently 
excellent fit completely fails to correctly predict the other 
two ratios measured for IC 10-SE. The observed modest 
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J=3-2 12 CO/ 13 CO ratio of eleven is four times lower than 
predicted, and the modelled J=4-3 12 CO intensity falls 
short of the observed value by a factor of two or more. 

Although this might suggest that single-component 
LVG analysis is physically irrelevant, this is not quite true. 
The kinetic temperatures and spatial densities implied by 
such fits often actually occur in the source. However, when 
they do, temperature and density generally do not refer to 
the same volume of gas, as will be clear from the following. 

4.1.2. Dual-component modelling 

For almost all of the sample galaxies, good fits based 
on two gas components could be obtained. Although the 
number of free parameters exceeds the number of inde- 
pendent measurements, this is a meaningful result because 
the actual combination of physically allowed parameters 
is sufficiently constrained. In order to reduce the number 
of free parameters, we have assumed identical CO isotopi- 
cal abundances for both gas components and assign the 
specific value [ 12 CO]/[ 13 CO] = 40. Different choices are 
possible, but reasonably small changes, for instance to ra- 
tios of 50 or 60 rarely lead to very different outcomes. 
We identified acceptable fits by searching a grid of model 
parameter combinations (covering the same range of tem- 
peratures, densities and gradients as described above) for 
matching model and observed line ratios, for various rela- 
tive contributions of the two components. 

We have rejected all solutions in which the denser gas 
component is also hotter than the more tenuous compo- 
nent, as we consider this physically unlikely on the large 
linear scales observed. From the remainder of solutions, 
we have selected characteristic examples and listed these 
in Tabled! The solutions are not unique, but delineate 
a range of values in particular parameter space regions. 
Variations in input parameters may compensate for one 
another, causing somewhat different input combination to 
yield identical line ratios, as indicated in the Table entries. 

The observed line ratios for the star-forming regions 
in all six galaxies listed in Table 5 require the presence 
of a rather dense component (typically njj 2 = 10 5 cm~ 3 ) 
but these components vary from fairly cool (IC 10-SE, 
NGC 6822-HV and probably NGC 5253) to warm (He2- 
10 and probably NGC 3077) to hot (NGC 1569). Column 
densities are well-established in the range A(CO)/dV = 
0.6-1.0 xl0 17 cm~ 2 /kms _1 . This dense component is as- 
sociated with a second component that is generally less 
dense, but hotter (typically 100 K). The actual density 
of this second component ranges from hh 2 = 100cm -3 
(IC 10-SE) to n H . 2 « 10 3 cm- 3 (NGC 1569, He 2-10 and 
probably NGC 3077). No signature of low-density gas 
is apparent in NGC 6822-HV and possibly NGC 5253. 
In these objects, the hot second component is still very 
dense (typically 10 5 cm~ 3 ). Note that the notable physi- 
cal difference between e.g. NGC 1569 and NGC 6822 is 
not immediately obvious as these objects have very sim- 
ilar 12 CO line ratios and identical J=l-0 12 CO/ 13 CO 



ratios. The important discriminators are the remaining 
12 CO/ 13 CO ratios. A relatively low J=3-2 12 CO/ 13 CO 
ratio (NGC 6822) is incompatible with a significant pres- 
ence of low-density gas at any temperature, whereas a 
high ,7=2-1 12 CO/ 13 CO (NGC 1569) requires the pres- 
ence of very hot and very dense gas, or just low-density 
gas, or a combination of both. Even then, we find it diffi- 
cult (as TableEl shows) to reproduce the very high J=2-l 
12 CO/ 13 CO apparently appropriate to NGC 1569. This 
suggests that the observations are in error, or that our as- 
sumptions are incorrect. Neither possibility can be ruled 
out. NGC 1569 is one of the most extreme (post) star- 
burst known (cf. Israel 1988) and its ISM appears to be 
subjected to intense processing (Lisenfeld et al. 2002). 
Thus, the less abundant 13 CO isotope may have selec- 
tively destroyed, so that its actual abundance is much less 
than assumed in our modelling. This would indeed lead to 
higher 12 CO/ 13 CO ratios. In the absence of further data 
we have chosen not to pursue this possibility, but it does 
suggest that (re) determination of the 13 CO intensities in 
NGC 1569 might be fruitful. 

4.2. Beam-averaged molecular gas properties 

The fraction of carbon locked up in the carbon monoxide 
molecule (i.e. the N(GO)/Nc ratio) as given by the chem- 
ical models by Van Dishoeck & Black (1988) is strongly 
dependent on the total carbon (and molecular hydrogen) 
column density. At column densities Nq > 10 18 cm -2 
practically all gas-phase carbon is in CO, whereas at col- 
umn densities Ac < 10 17 cm -2 essentially all carbon is in 
atomic form. Using these chemical models and the radia- 
tive transfer model parameters summarized in Tabled we 
have estimated the beam-averaged carbon monoxide col- 
umn densities as well as the beam-averaged column den- 
sities of all (atomic and molecular) carbon present in the 
gas-phase. Total carbon column densities Ac may be con- 
verted into total hydrogen column densities Ah as a func- 
tion of the [C]/[H] abundance and the fraction of carbon in 
the gas-phase. For galaxies with metallicities below solar, 
the [C]/[H] abundance can be estimated from the [C]/[0] 
versus [0]/[H] diagrams given by Garnett et al (1999) and 
the [0]/[H] metallicities summarized in Table0J We fur- 
thermore assume a gas-phase carbon-depletion factor 5c 
= 0.27. By subtracting observed neutral hydrogen column 
densities A (HI) from the total hydrogen column densities 
thus obtained, beam-averaged molecular hydrogen column 
densities are determined. Finally, these are used to ob- 
tain total molecular hydrogen masses and estimates of the 
CO-to-H2 conversion factor X. Tabled gives the results 
for those galaxies where this procedure was meaningful. 
We find that at least in the low-metallicity starburst ob- 
jects, overall CO column densities are low (which is re- 
flected by the relatively low velocity-integrated intensi- 
ties in Table|2J). We also find that in these objects, only a 
small fraction of all carbon is in CO and that most carbon 
should be in neutral or ionized atomic form. As a conse- 
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Table 7. Beam-averaged physical parameters 



Galaxy 


Abundances 
N c /N (CO) Nh/Nc 


Column Density 
JV(HI) iV(H 2 ) 
(cm -2 ) (cm~ 2 ) 


Mass 
M(H 2 ) 
(Ma) 


CO-H 2 Ratio 
X 

(cm" 2 /Kkms _1 ) 


Haro 2 


21±2 


4 x 10 4 


1.6 x 10 21 


5 x 10 21 


7 x 10 7 


9 x 10 20 


IC 10-SE 


5±1 


1 x 10 5 


3.6 x 10 21 


24 x 10 21 


3 x 10 6 


12 x 10 20 


NGC 1569 


46±23 


1 x 10 5 


3.7 x 10 21 


11 x 10 21 


1 x 10 7 


72 x 10 20 


He 2-10 


11±2 




1.9 x 10 21 






< 4 x 10 20 


NGC 6822-HV 


22±3 


6 x 10 4 


1.6 x 10 21 


10 x 10 21 


5 x 10 5 


60 x 10 20 



Note: HI column densities from Bravo- Alfaro et al. (2004: Haro 2), Madden et al. (1997: IC 10), Stil & Israel (2002: NGC 1569), 
Kobulnicky et al. (1995: He 2-10), Israel et al. (2003: NGC 6822). 



Table 8. Environment and molecular gas parameters 



Galaxy 


(2-i)/(i-o) 


12 co 

(3-2)/(2-l) 


12 GO/ 13 

(3-2)/(l-0) (1-0) (2-1) 


CO 
(3-2) 


N c /N(GO) 


log N(CO) 
(cm- 2 ) 


log AT(H 2 ) 
(cm" 2 ) 


Compact Galaxies 


Type I (Starburst) 
Type II (Active) 
Type III (Quiet) 


1.1 

0.9 
0.5 


1.1 
0.8 


1.2 20 23: 
0.7 17: 


13 


21 
21 


16.6 
16.2 


22.2 
21.7 


Spiral Galaxy Centers 


Starburst 
Quiet Nucleus 


1.0 
0.5 


0.7 
0.6 


0.7 10 10 
0.3 7 6 


11 


3 
6 


18.1 
18.0 


21.7 
21.6 


LMC/SMC Clouds 


Star-Forming 
Quiet 


1.1 
0.7 


0.9 
0.5 


1.0 12 8 
0.4 10 14 











qucncc, molecular hydrogen densities are not proportion- 
ally lower. Low integrated CO intensities and relatively 
normal molecular hydrogen column densities together im- 
ply the high X-values given in TableEl 

4.3. Molecular gas in compact galaxies 

Although the database is admittedly still limited, we can 
make a few statements as to the general nature of molecu- 
lar gas in compact and dwarf irregular galaxies. On the ba- 
sis of CO spectroscopy we may distinguish different types 
of environment in compact galaxies, with properties sum- 
marized in Table|S] We have included results obtained in 
a similar way for a sample of high-metallicity starburst 
galaxy centers (NGC 253, IC 342, Maffei 2, M 83, and 
NGC 6946 - Israel et al. 1995, Israel & Baas 2001, 2003) 
and for a quiescent centers such as occur in NGC 7331 
(Israel & Baas 1999) and the Milky Way (Bennett et 
al. 1994). The former are characterized by CO intensi- 
ties slowly decreasing with increasing rotational transi- 
tion, as well as 12 CO/ 13 CO ratios of about 10. The lat- 
ter are weaker CO emitters, more difficult to detect and 



less frequently included in surveys. They have CO inten- 
sities much more rapidly decreasing with rotational tran- 
sition, and they have lower 12 CO/ 13 CO ratios around 6 
- 7. Finally, we also included average ratios obtained for 
both quiescent and actively star-forming molecular clouds 
in the low-metallicity environments of the Large and the 
Small Magellanic Clouds (Bolatto et al. 2000, 2005). 

Type I (starburst compact) has CO emission with 
J T m bdl^ rising with increasing J-transition: both the (2- 
1)/ 1—0) and the (3-2)/(2-l) ratios are close to or even 
exceeding unity. At the same time, the 12 CO/ 13 CO iso- 
topical ratios are quite high, typically of the order of 20. 
The CO results obtained towards NGC 1569, He 2-10, 
NGC 5253 and NGC 6822-HV are all in this category, as 
are those of the object UM 465 mentioned earlier. The 
molecular clouds sampled are dominated by intense star- 
formation activity. This may be representative of a large 
part of the galaxy (NGC 1569, He 2-10) or only of a spe- 
cific location (NGC 6822-HV). The rotational line ratios 
imply high overall temperatures, and the isotopical ratios 
are characteristic of relatively low molecular gas beam- 
averaged column densities. This is clearly illustrated in 
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TablcOD The line ratios of Type I compact galaxies are 
very similar to those of spiral galaxy starburst centers 
and the low-metallicity starforming clouds in LMC and 
SMC, except for much higher 12 CO/ 13 CO ratios. This is 
characteristic of similar (high) temperatures but lower op- 
tical depths in the starburst dwarf galaxies. Hence, in the 
compact galaxies beam-averaged CO column densities are 
more than an order of magnitude lower than those in star- 
burst galaxy centers and star- forming LMC/SMC clouds. 
In the compact galaxies most carbon is in atomic form. 
The beam-averaged molecular hydrogen columns are very 
similar to those of the starburst centers. The line ratios 
of IC 10-SE are different from any other set. They most 
resemble a mixture of active and cold LMC/SMC clouds. 

Because integrated CO intensities are low and the 
H 2 column densities implied by the analysis are not, 
Y-factors in compact galaxies very significantly exceed 
those of starburst centers. We have calculated values X = 
(12,72,60) xlO^cm-^Kkms- 1 )- 1 for IC 10, NGC 1569 
and NGC 6822 respectively. These values are in good 
agreement with conclusions reached earlier and indepen- 
dently for IC 10 by Madden et al. (1997) and for NGC 1569 
by Lisenfeld et al. (2002). They also lend further support 
to the strong metallicity dependence of the Y-factor pro- 
posed by Israel (1997, 2000). 

Type II (active compact) has Haro 2 as a prototype. 
The line parameters are similar to those of Type I ob- 
jects, but the rotational line ratios are lower and less than 
unity. The difference is clearly seen in the 12 CO (3-2)/(l- 
0) ratio which is almost half that of the Type I average. 
Most of the galaxies in the top panel of Tabled] includ- 
ing NGC 3353 (Haro 3) are of Type II. They are actively 
forming stars, but the CO spectroscopy is not dominated 
by star-forming regions as in the Type I objects. IC 10 
and NGC 6822, now included in Type I, would probably 
fall in this category if we had sampled the whole galaxy 
instead of only the most actively star- forming part. 

Even for Haro 2 the observational data are insufficient 
to allow a meaningful two-component analysis, so that the 
derived parameters are based only on an overly simplified 
single-component analysis. Nevertheless, as in Type I ob- 
jects, carbon monoxide appears to be only w 5% of all 
gas-phase carbon. Carbon monoxide columns appear two 
orders of magnitude lower than those in starburst nuclei, 
but molecular hydrogen columns are virtually identical. 
For Haro 2 we find X = 9x 10 20 cm _2 (Kkms -1 ) -1 practi- 
cally identical to the value derived in a different manner by 
Barone et al. (2000). Although this is much lower than the 
average derived for the Type I objects, it is very close to 
the value expected for this metallicity in the diagrams by 
Israel (1997, 2000). If we assume, just for the sake of argu- 
ment, the same metallicity for all type II galaxies, we can 
perform an analysis similar to that of Haro 2 also for the 
other galaxies in this particular subsample. Not surpris- 
ingly, we obtain very similar CO, C and H 2 columns (with 
a mean log Y(H 2 ) = 21.8). For NGC 2537 and NGC 2976 
we would find X = 10 - 30 x 10 20 cm~ 2 / K kms" 1 , but for 
the much more distant galaxies NGC 4194 and NGC 6052 



we would obtain X values an order of magnitude lower, 
not substantially different from the Solar Neighborhood 
value X = 2 x 10 20 cm" 2 /Kkms" 1 . The fact that these 
are the most luminous galaxies in the sample (cf. TablcEJ 
is probably more relevant to this result than the fact that 
they are the most distant galaxies. 

In the Type II objects, the 12 CO (3-2)/(l-0) ratio is 
practically constant (at about 0.67) over the full range 
of infrared colors / 60 //ioo- This mean 12 CO (3-2)/(l-0) 
ratio is, in fact, identical to the mean of 0.66 found by 
Yao et et al. (2003) for a sample of 60 infrared-luminous 
galaxies (with individual values varying from from 0.22 to 
1.72). Yao et al. also presented evidence suggesting X to 
be lower than the Solar Neighbourhood value by an order 
of magnitude, making their sample more similar to our 
starburst centers than to our compact galaxies. 

Type III (quiet compact) consists of those galaxies 
where low (2-1) /(1-0) line ratios ~ 0.5 suggest that on 
the whole little is going on. Their molecular gas is not 
very dense, and may also not be very warm. 2 Zw 40 is 
a member of this category, weak in J = 1 — 0, barely de- 
tected in J = 2 — 1, and not in J = 3 — 2. This class 
also apears to include NGC 2730, NGC 3659, NGC 4532, 
NGC 6570, NGC 7732, and IC 3521 from the sample ob- 
served by Albrecht et al. (2004). It would be of interest 
to determine further line ratios in order to constrain the 
physical conditions of the molecular gas, but the weakness 
of emission in all but the J = 1-0 transition makes this a 
difficult task. 

5. Conclusions 

1. We have observed 12 CO line emission in the J=2— 1 
and .7=3-2 transitions from 1 1 compact (dwarf) galax- 
ies. In four cases, limited maps were made. 

2. By combining our data with those from the literature, 
we eliminated beam-dilution effects and established ac- 
curate line ratios for the first three 12 CO rotational 
transitions in a sample of a dozen objects, and lim- 
ited the range of possible values in another four. For 
six objects, 12 CO/ 13 CO isotopical ratios were culled 
from the literature. 

3. Radiative-transfer (LVG) modelling shows that in 
most of the observed galaxies warm molecular gas oc- 
curs with temperatures typically = 50 — 150 K as 
well as dense gas with n(H 2 ) > 3000 cm~ 3 . Models 
using two distinct molecular gas components produce 
results different from and better than those obtained 
by modelling only a single component. 

4. Chemical modelling implies that in the observed galax- 
ies only a small fraction (typically 5%) of all gas-phase 
carbon is in the form of CO, the remainder being in 
the form of neutral carbon ([CI]) and especially ionized 
carbon ([CII]). As a consequence, CO column densities 
are quite low and of the order of 10 16 cm~ 2 . 

5. Molecular hydrogen column densities are high, of the 
order of 10 22 cm~ 2 , and confirm the large CO-to- 
H 2 conversion factors, in the range X = 10 21 — 
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10 22 cm~ 2 / K kms -1 , found earlier for low-metallicity 
environments by different methods. 
6. The CO spectroscopy of compact galaxies may be clas- 
sified into three different types. Type I is characterized 
by high rotational and isotopic ratios, and reflects hot 
and dense molecular clouds dominated by star-forming 
regions. Type II has lower ratios, in particular the 
12 CO (3-2)/(l-0) ratio is much lower than in Type I, 
and identical to the mean found for infrared-luminous 
galaxies in general. Type III has a low 12 CO (2-1)/ (1- 
0) ratio indicative of not very dense and possibly rel- 
atively cool molecular gas, and appears to represent 
quiescent compact galaxies. 
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